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ABSTRACT 

Context. Water is present during all stages of star formation: as ice in the cold outer parts of protostellar envelopes and dense inner 
regions of circumstellar disks, and as gas in the envelopes close to the protostars, in the upper layers of circumstellar disks and in 
regions of powerful outflows and shocks. Because of its key importance in the understanding of its origin in our own solarsystem, 
following the evolution of water all the way to the planet-forming disk is a fundamental task in research in star formation and astro- 
chemistry. 

Aims. In this paper we probe the mechanism regulating the warm gas-phase water abundance in the innermost hundred AU of deeply 
embedded (Class 0) low-mass protostars, and investigate its chemical relationship to other molecular species during these stages. 
Methods. Millimeter wavelength thermal emission from the para-Hj^O 3i,3 - 22.0 (-£ii=203.7 K) line is imaged at high angular reso- 
lution (0'.'75; 190 AU) with the IRAM Plateau de Bure Interferometer towards the deeply embedded low-mass protostars NGC 1333- 
IRAS2A and NGC 1333-IRAS4A. 

Results. Compact HJ,^0 emission is detected towards IRAS2A and one of the components in the IRAS4A binary; in addition 
CH3OCH3, C2H5CN, and SO2 are detected. Extended water emission is seen towards IRAS2A, possibly associated with the out- 
flow. 

Conclusions. The results complement a previous detection of the same transition towards NGC 1333-IRAS4B. The detections in 
all systems suggests that the presence of water on < 100 AU scales is a common phenomenon in embedded protostars and that the 
non-detections of hot water with Spitzer towards the two systems studied in this paper are likely due to geometry and high extinction 
at mid-infrared wavelengths. We present a scenario in which the origin of the emission from warm water is in a flattened disk-like 
structure dominated by inward motions rather than rotation. The gas-phase water abundance varies between the sources, but is gener- 
ally much lower than a canonical abundance of 10""*, suggesting that most water (>96%) is frozen out on dust grains at these scales. 
The derived abundances of CH3OCH3 and SO2 relative to Hj^O are comparable for all sources pointing towards similar chemical 
processes at work. In contrast, the C2H5CN abundance relative to Hj^O is significantly lower in IRAS2A, which could be due to 
different chemistry in the sources. 

Key words, astrochemistry - stars: formation - planetary systems: protoplanetary disks - ISM: abundances - ISM: individual (NGC 
1333-IRAS2A, NGC 1333-IRAS4A) 

1 . Introduction Bure Interferometer (PdBI), one of the few water lines that can 

,., . . , , • • ,- r • T ,. be routinely observed from Earth. The aim of the observations 

Water is a key molecule m regioiis of star formation. In cold ^^ ^^ determine the origin of water emission lines in the inner 

5_] outer regions of protostellar envelopes and dense mid-planes loO AU of these protostars. 

C^ of circumstellar disks, it is the most abundant molecule in the t,. j.j ujj.» j 

. . . , . , iTTTT-; n u^^j it: n u^^o t The young, dust- and gas-enshrouded protostars provide an 

ices on dust grains (e.g. iGibb et al.ll2004l: iBoogert et al.ll2008l) . . .iu . f , a- f .u . f .■ 

- ?' . % ' " — ; — " — 7-* " — rr important laboratory for studies of the star formation process. 

In regions of increased temperatures (e.g., due to protostellar ^ ,, . j .u j .u v 1 u . .u . n 

, .'^ . ,,N., ,,, ■ On the one hand, they provide the link between the prestellar 

heating or active shocks) it desorbs and becomes an impor- , ^, -^ J'^ , . ^ ,vcr\ \ a ^u ^u 

° , J, r=; —^. Tmz^Ti w^ — r 1 1 n -.r — 7^ — I cores and the young stellai" objects (YSOs) and on the other 

tant coolant for gas (ICeccarelh et al.lll996t INeufeld & KaufmanI ■ j ■» ■ n 1 »u » j ■ »u ^ »u ■ ■*• 1 j»- e 

nrrr-ri . ° ^ ^ : " • ' : — : — \ . _ ^ hand it is likely that during these stages the initial conditions for 

I1993D . As one of the major oxygen reservoirs it also innu- .. 1 . . u • 1 1 .• ?.u * » n 1 

, , , , • ^ I T^T-; , n „, , I the later astrochemical evolution of the protostellar envelopes 

ences the oxygen based chemistry (e.g. Ivan Dishoeck & Blakd j j- 1 » ui- u j / \r- » 1 oaaoa t »u i- 

i --.__ i ^, . -"^ , • , ^ \ — '■ ; and disks are established (e.g., Visseret al. 2009). In the earli- 

119981) . This paper presents high-angular resolution observa- , , » ^ .u a a ^ ■ .1 1 

T^ — rr^iir^ 5,„o ^ ^,, j- , ,,is^ • , , cst stagcs most of the gas and dust is present on large scales 

tions (C.'TS) at 203.4 GHz of the H'**0 isotopologue towards • ,. . . n 1 u .u . .1 

,,,,,,,_, ^^ , 2 fa jjj (jjg protostellar envelopes where the temperatures are low 

two deeply embedded (Class 0) low-mass protostars using the ,,n onV\ a » 11 ■ ^ a- ^ iij 

, ^ ,. •»/•„• , ■ /T„ . , ,^ ™ , (10-20 K) and most molecules, including water, are locked up 
Institute de Radioastronomie Millimetnque (IRAM) Plateau de . ■ m »u 1 \. a ^. ^- c 
-1 \ y ^g j(.g Qjj grains. Nevertheless, space-based observations from 

* Based on observations carried out with the IRAM Plateau de Bure ISO, SWAS, Odin, Spitzer and Herschel show that significant 
Interferometer. IRAM is supported by INSU/CNBRS(France), MPG amounts of gas-phase water are present towards even deeply em- 
(Germany) and IGN (Spain). bedded protostars, pointing to significant heating and processing 
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of water-ice covered grains (iNisini et aL 1999t Ceccarelli et alJ 



1999"; 'Maret et al."2002'; 'Bergin et al."2003' Franklin et a l.ll2008 
Watson et al. 2007; van Dishoeck et al. 2011; Kristense n et alj 



20101) . 



Exactly which mechanism regulates where and how the heat- 
ing and evaporation of the water ice takes place is still debated. 
Proposed locations and mechanisms include the innermost enve- 
lope, heated to temperatures larger than 100 K by the radiation 
from t he central protosta r on scales <100 AU (ICeccarelli et al.l 
1 1998'; Maret et al. 7002), in shocks of the protostellar outflow 
(Nisinietal. 1999), or in the disk - heated either by shocks 
related to ongoing accretion onto the circu mstellar disks or in 
their surface layer s heated b y the protostar (IWatson et al.ll2007t 
lJ0rgensen & van D ishoeck 2010al) . 

Recently, Watson et al. (2007) observed 30 Class proto- 
stars, including the sources in this paper, with the infrared spec- 
trograph (IRS) on board Spitzer. The protostar NGC 1333- 
IRAS4B (IRAS4B) showed a wealth of highly excited H2O 
lines, and was the only source in the sample in which the wa- 
ter lines were detected. By modeling the spectra they concluded 
that because of the high inferred density and temperature of the 
water, it originates in the vicinity of the protostar, likely accre- 
tion shocks in the circumstellar disk. However, the Spitzer data, 
and similar observations from space, lack the spatial and spec- 
tral resolution to unambiguously associate the emission with 
the components in the protostellar systems - and can for the 
shorter wavelengths also be hampered by extinction by the dus t 
in the protostellar envelopes and disks. Herczeg et al.l (l2011h 
analyzed the emission from warm water lines in maps toward 
IRAS4B obtained with the PACS instrument on the Herschel 
Space Observatory. The maps show that the bulk of the warm 
water emission has its origin in a outflow-driven shock offset by 
about 5" (1200 AU) from the central protostar. An excitation 
analysis of the emission furthermore shows that the temperature 
in the shocked region is about 1500 K, which can simultane- 
ously explain both the far-IR (Herschel) and mid-IR (Spitzer) 
water lines. 

Higher angular and spectral resolution is possible through 
interferometric observations of Hj^O which has a few 
transitions that can be imaged from t he ground (e.g. 



JacQ et al. 1988; Gensheimer et al. 1996t Ivan der Tak et al.l 
200 6). J0rge nsen & van Dishoeck (2010ai) report observations 
of IRAS4B' with die PdBI. They detect the H^^O 3 1,3 - 22,0 hne 
at 203.4 GHz, and also resolve a tentative velocity gradient in 
the very narrow line (~2 km s"' width at zero intensity), indicat- 
ing an origin of the Hj^^O emission in a rotational structure. The 
H2O column density is around 25 time s higher than previo usly 
estimated from mid-IR observations bv lWatson et al.l(l2007l) . on 
similar size scales (inner 25 AU of the protostar). 

This study presents similar high angular resolution Hj^O 
millimeter wavelength observations for the sources NGC 1333- 
IRAS2A (IRAS2A) and NGC 1333-IRAS4A (IRAS4A), which 
are compared to IRAS4B. The three sources are all deeply em- 
bedded Class protostars in the nearby NGC 1333 region of the 
Perseus molecular cloucfl identified through earl y infrared and 
(sub) r ninimeter waveleng th observations (e.g., iJennings et alj 
11987"; 'Sandell et al."l99l'). All sources were mapped by Spitzer 
(J0rgensen et al. 2006; Gutermuth et al. 2008), and were part 
of the high resolution (~1") PROS AC surv ey of Class O/I 
sources using the Submillimeter Array (SMA; |j0rgensenet alJ 



l2007i 120091) . Their similar luminosities (20, 5.80, 3.8 Lq for 
IRAS2A, IRAS4A and IRAS4B, respectively) and envelope 
masses (1.0, 4.5^, 2.9 Mq, respectively) from those surveys, 
and thus presumably also similar evolutionary stages, make them 
ideal for studying and comparing deeply embedded young stellar 
objects in the same environment. 

IRAS4A is a binary system with a 178 separation 
(Looney et al. 2000; Reipuith et al. 2002; Girart et al. 2006). as- 



socia t ed with a larger scal e jet and outflow (e.g. iBlake et al.l 
Il995l : |j0rgensen et al.l2007b . The jet axis has been deduced from 
water maser observations associated to IRAS4A-NW to 2° in 
the plane of the sky and at a position angle (PA) of -50 ° (north 
through east. Marvel et al. 2008; Desmurs et aUl20p6|). almost 
perpendicular to the CO outfl ow axis (PA= 0°-45°; lBlake et all 
I1995I; lj0rgensen et ani2007h . The mas s of the NW sour ce is es- 
timated from modeling to 0.08 Mq (IChoi et al .11201 Oh . no es- 
timate of the mass of the SE source exists. The envelope is 
dominated by inward motion, as evident from the inverse P- 
Cygni profiles in the spectra and modeling of gas kinematics 
(Y)\ Francesco et all 1200 it lAttard et al.l 120091; iKristensen eTal] 
201(J 

The other source, IRAS2A, shows two bipolar outflows , 
one sheU-Hke directed N-S (PAsj -h25°, lj0rgensen etd]l2007l) 
and the other a jet-like directed E-W (PA^ -Hl05 °,| Liseau et all 
Il988t ISandeUetal.ll 1994 lj0rgensenetani2OO4l) . Water maser 
emission has been detected towards I RAS2A, partly asso ciated 
with the northern redshifted outflow (iFuruva et al.l 120031) . One 
proposed solution to the presence of a quadrupolar outflow is 
that IRAS2A is a bin ary system with 0'.'3 (75 AU) separation 
(lJ0rgensen et al.ll2004l) . 

The paper is laid out as follows. Sect. |2]presents the interfer- 
ometric observations of the three Class sources that the paper 
is based upon. Sect. [3] presents the results from the analysis of 
the spectra and maps. Sect. |4| discusses the detection of water, 
and the column density estimates. Sect.|5]presents a summary of 
the important conclusions and suggests possible future work. 



2. Observations 

Three low-mass protostars, IRAS2A, IRAS4A-SE and -NW in 
the embedded cluster NGC 1333 in the Perseus molecular cloud 
were observed using the PdBI. IRAS2A and IRAS4A were ob- 
served in two configurations each: the B configuration on 6, 15 
and 17 March 2010 and the C configuration on 9 December 
2009. The sources were observed in track-sharing mode for a 
total of 12 hours in B configuration and about 8 hours in C con- 
figuration. Combining the observations in the different config- 
urations, the data cover baselines from 15.4 to 452 m (10.5 to 
307 k/l). The receivers were tuned to the para-Hj^O 3 1,3 - 22,0 
transition at 203.40752 GHz (1.47 mm) and the correlators were 
set up with one unit with a bandwidth of 40 MHz (53 km s"') 
centered on this frequency providing a spectral resolution on 460 
channels of 0.087 MHz (0.115 km s') width. 

The data were calibrated and imaged using the CLIC and 
MAPPING packages, which are parts of the IRAM GILDAS 
software. The general calibration procedure was followed: reg- 
ular observations of the nearby, strong quasar 0333-H3221 
were used to calibrate the complex gains and bandpass, while 
MWC349, 0333H-3221 and 3C84 were observed to calibrate the 
absolute flux scale. Integrations with significantly deviating am- 
plitudes/phases were flagged and the continuum was subtracted 



In t his paper we adop t a distance of 250 pc to the Perseus molecular 



In t nis paper we adop t a distance ot /du pc to 
cloud l lEnoch et al.ll2006| . and references therein) 



The luminosity and mass estimates for IRAS4A refer to the binary 
system. 
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before Fourier transformation of the data. The resulting beam 
sizes using natural weighting, along with other data, are given 
in TableB The field of view is 25"(HPBW) at 1.45 mm and the 
continuum sensitivity is limited by the dynamical range of the 
interferometer while the line data RMS is given in Table [1] The 
uncertainty in fluxes is dominated by the flux calibration error, 
typically 20%. 



Table 1. Observation parameters. 



Source 



Beam" 



PA 



RMS" 



IRAS2A 
IRAS4A(-SE/NW) 



a'87xa'72 

a'86xa'7o 



63.5° 
63.4° 



16.1 
16.5 



Major X minor axis. ' ' Units of mjy beam channel 



3. Results 

Continuum emission from the targeted sources (IRAS2A and the 
IRAS4A-NW and IRAS4A-SE binary) was detected at levels 
expected from previous high resolution millimeter wavelength 
continuum observations: parameters from elliptical Gaussian fits 
to the continuum emission for the sources are given in Table |2] 
The continuum emission toward s the sources is resolved and the 
fluxes concur with the results of |j0rgensen et al.l (|2007|) . assum- 
ing a power-law spectrum (Fy oc V) with a~2.3-2.7, expected 
for thermal dust continuum emission. 

Fig. [T] shows the spectrum towards the continuum peak of 
these s ources, as well as IRAS4B (J0rgensen & van Dishoecy 
1201 Oah . A number of lines are detected towards IRAS2A and 
one component in the IRAS4A protobinary system, IRAS4A- 
NW. IRAS4A-SE in contrast does not show any emission lines 
in the high-resolution spectra, and it is therefore only listed in 
the table showing the characteristics of the continuum emission . 

The Jet Propulsion Laboratory (JPL, iPickett et all 1 19981) 
and the Cologne Database for Molecular Spectroscopy 
(CDMS. lMuller et aP UoOl) databases were queried through the 
Splatalogue interface to identify the lines. The line widths and 
strengths vary between the sources and because of the broader 
line widths some blending between the different CH3OCH3 lines 
occurs in both 1RAS2A and 1RAS4A-NW. As seen in Fig.fT] sev- 
eral lines are detected: the target water line, SO2 and several lines 
of CH3OCH3. C2H5CN is seen in IRAS 4A-NW and 1RAS4B, 
but not in 1RAS2A. 

Particular interesting for this study is the assignment of 
the H^'^O 3 1,3 - 22,0 transition (£'u=203.7 K). No other fines 
faU witfiin ±1 km s"' of tfie Hj^O 3 1,3 - 22,0 line listed 
in tfie spectral line catalogs. Tfie lines tfiat fall close (+1- 
2 km s-i) are 13CH2CHCN, 13CH3CH2CN and (CH3)2CO. Botfi 
tfie '^CH2CHCN and '^CH3CH2CN isotopologues fiave otfier 
transitions wfiicfi would fiave been detected, wfiereas a possi- 
ble transition of (CH3)2CO fias a too low intrinsic line strengtfi 
to be plausible. 

In Fig. 12] tfie integrated intensity (zerotfi order moment) and 
velocity (first order moment) maps for tfie Hj^^O transition are 
sfiown for tfie observed sources. For 1RAS2A tfie calculations 
were made over tfie velocity interval from 4.4 to 1 1 .4 km s"' and 
for 1RAS4A-NW from 4.7 to 9.6 km s '. Tfie integrated inten- 
sity maps for 1RAS2A sfiow some extended emission, wfiereas 
for 1RAS4A-NW tfie emission is compact around tfie continuum 
peak witfi a sligfit elongation in tfie nortfi-soutfi direction. Tfie 
velocity maps for botfi 1RAS2A and 1RAS4A (rigfit fiand panels 
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Fig. 1. Spectra extracted in a pixel towards tfie continuum peak 
positions for 1RAS2A, 1RAS4A-NW and 1RAS4B. Tfie identi- 
fied lines are marked, including tfie targeted H^^O 3 1,3 -22,0 line. 
Parameters from tfie fits are tabulated in Table IA]2. Tfie numbers 
above tfie lines and tfie colors correspond to tfie numbers (emis- 
sion lines) and colors in Fig.|3] Table[3] Table lA.fl and TablelAl2. 
Tfie spectra fiave been binned to a widtfi of 0.23 km s"' (i.e., 2 
cfiannels per bin), and Gaussian profiles fiave been fitted to tfie 
identified lines. Tfie systemic velocity of tfie H2'^0 line is marked 
by tfie vertical dotted line at visi = 7.0 km s"'. Note tfie differ- 
ence in scale of tfie y-axis for tfie spectra. 



in Fig.|2| sfiow no obvious spatial gradients in velocity, neitfier 
on source nor along tfie outflow in 1RAS2A. 

For estimates of tfie pfiysical conditions in tfie gas, circular 
Gaussian profiles were fitted to tfie emission of eacfi line in tfie 
(m, v)-plane. Tfie fits were performed to tfie line intensities typi- 
cally integrated over tfie line widtfis (v'lsi+FWHM) and are given 
in Table [3] Table IaTI lists tfie molecular parameters adopted in 
tfiis paper of tfie target line Hj^O 3 13 - 22,0, and otfier detected 
lines, wfiile Table lAl2 lists tfie position, FWHM and integration 
interval for tfie different lines. 



4. Analysis and discussion 

Tfie detected water lines provide important constraints on tfie 
pfiysical and cfiemical structure of tfie deeply embedded low- 
mass protostars, wfiicfi we discuss in tfiis section. Tfie first part 
of tfie discussion (Sect. 14.1114.31 ) focuses on tfie detection of 
warm water on scales comparable to tfie protostellar disk and tfie 
extended emission found in tfie outflow towards 1RAS2A. Tfie 
implications for tfie geometry and dynamics are also discussed. 
Tfie second part (Sect. l4.4H4.6b discusses excitation conditions 
in tfie water emitting gas and tfie column densities and relative 
abundances of water and otfier species detected on small scales. 



M. V. Persson et al.: Subarcsecond resolution observations of warm water towards three deeply embedded low-mass protostars 

Table 2. Parameters for the sources from elliptical Gaussian fits to their continuum emis- 
sion. EiTors given are the statistical eiTors from the fitting routine. 



IRAS2A 



IRAS4A-NW 



IRAS4A-SE 



Flux (±)'' 0.24 (0.002) Jy 

R.A i±)'' 03:28:55.57 (0'.'004) 

Dec (±)* 31:14:37.03 (ft' 004) 

Extend I'.'lx l'.'0(23°) 



0.88 (0.025) Jy 
03:29:10.44 (a'02) 
31:13:32.16 (a'02) 
2'.'0 X 1'.'2 (-53°) 



1.46 (0.021) Jy 
03:29:10.53 (a'005) 
31:13:30.92 (a'006) 
l'.'4x l'.'0(50°) 



'"' Statistical error from fit, the calibration dominates the uncertainty in the flux determination. 
**' J2000 epoch. Statistical error from fit. *^' Major x Minor axis and position angle for elliptical 
Gaussian profile in (ii, v) plane. The statistical error of the extent is typically 0.5 - 3%. 
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Fig. 2. Integrated intensity (left) and veloc- 
ity (right) maps for the H^'^O 3 13 - 22,0 line 
in IRAS2A (from 4.4 to 'l 1.4 km s"') and 
IRAS4A (from 4.7 to 9.6 km s"'). Contours 
show integrated intensity and are in steps of 4cr 
for the intensity map, 8cr for the velocity map, 
and start at 3cr. The position of the continuum 
emission is marked by a red plus sign, an out- 
flow by the solid line and maser emission by 
the dashed line. Positions where spectra have 
been extracted in IRAS2A along the outflow 
are marked with letters a and b and crosses with 
magenta and green colors in the integrated in- 
tensity map. 



4. 1 . Water on small scales 

The key result from the presented observations is the detection 
of compact H^^O in IRAS2A and I RAS4A-NW , addin g to th e 
discovery towards IRAS4B ( Jdrgens en & van Dishoeckii2010ah . 
Together, the observations demonstrate the presence of water 
emission on 100 AU scales - much sma ller than probe d by the 
recent Herschel (iKristensen et al.l [ZOlOh and ISO (Nisi ni et aP 
11999 ) observations. 

It is particularly noteworthy that the water emission is 
present in all these systems - and the integrated emission (spa- 
tially and spectrally) is in fact stronger in both IRAS4A and 
IRAS2A than in IRAS4B. This is in contrast to the space ob- 
servations where the Herschel HIFI data show equally strong 
or stronger e mission lines in IRAS4B than the other sources 
(IKristensen et al. 201 0) and the Spit zer data only show detec- 
tions in IRAS4B ( Watson etalJI 20071) . 

The extent of the emission estimated from the sizes of the 
circular Gaussian fits in the (m, v) plane (TableO differ between 
the sources: for IRAS2A emission from the different species is 
generally more extended (O'.'VO and 1'.'45; 87-181 AU radii) than 
in IRAS4A-NW (a.'37 and 0'.'88; 46-110 AU radii). Still, for 
both sources (as well as IRAS4B), these results show that the 
bulk of the emission is truly coming from the inner 50-100 AU 
of the protostellar systems. 



The positions of the emission differ slightly: Fig.|3]shows the 
peak positions of Gaussian fits in the (M.vj-plane to the integrated 
intensity of the detected lines relative to the continuum position 
for the two sources, IRAS2A and IRAS4A-NW, from this paper 
as well as IRAS4B. The size of the ellipses show the statistical 
uncertainty in the Gaussian fits, and the numbers correspond to 
those in Fig. [J] Table [3] Table [Al] and Table |A]2 identifying 
each line. The position of the lines in the different sources show 
a systematic offset, in IRAS4A-NW and IRAS4B the lines are 
all (except C2H5CN in IRAS4B) offset about 0'.'2 west of the 
continuum peak. When fitting the weak line of CH3OCH3 (line 
6) in IRAS2A, an extra Gaussian was added, centered on the 
stronger extended emission to minimize its affect on the fit of the 
line. The offsets are smaller than the beam, but still significant 
given the signal-to-noise of the data. That the shifts are physical 
is further strengthened by the observation that they are similar 
for the observed species toward each source. 



IChandler et al.l (l2005h detect similar offsets in the line emis- 
sion of complex organic molecules relative to the continuum 
source towards IRAS 16293-2422 combining high-resolution 
SMA observations with archival VLA data. They argue that an 
outflow shock eastward of one component in the protostellar bi- 
nary ("source A") is responsible for enhancing the abundances, 
and sputtering the ices off the grains. Because the offsets seen 
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Table 3. Parameters for the lines detected in the sources from circular Gaussian fits to the integrated 
line maps. Molecule specific data are tabulated in Table lA.fl and the integration interval (in km s ') in 
Table lAl2 in the appendix. 



No. Molecule Transition 



Flux" 
(Jy km s"') 



IRAS2A 
Offset" 

(") 



Size^ 

(") 



Flux" 
(Jykm s"') 



IRAS4A-NW 
Offset" 

(") 



Size^ 
(") 



1 


CH3OCH3 


33,0 " 22.1 EE 


0.47" 


-0.02; -0.002 


0.70 ± 0.09 


0.44" 


-0.21 


-0.04 


0.70 ± 0.09 


2 


CH,OCHi 


33.0 - 22.1 AA 


_d 


-;- 


- 


_rf 


- 


- 


- 


3 


CH,OCHi 


33.0-22,1 AE 


0.8P 


-0.16; -0.32 


1.45 ±0.10 


0.26'' 


-0.14 


-0.03 


0.57 ±0.12 


4 


CH3OCH3 


33.1-22, EE 


-'■ 


-;- 


- 


_e 


- 


- 


- 


5 


Hi^O 


3 1.3 -22,0 


0.98 


-0.07; -0.12 


0.83 ± 0.05 


0.27 


-0.18 


-0.07 


0.61 ±0.12 


6 


CH3OCH3 


33.1- 22,1 EA 


0.30 


-0.35; -0.52 


1.23 ±0.22 


0.06 


-0.25 


-0.10 


0.44 ± 0.26 


7 


CHsCN 


232,22 - 222,21 


< 0.04^ 


-; - 


- 


0.27 


-0.20 


-0.05 


0.37 ±0.11 


8 


SO2 


12o,12 - 111, 11 


2.73 


-0.07; +0.003 


1.00 ±0.02 


1.22 


-0.22 


-0.01 


0.88 ± 0.05 



*"' Flux from fitting circular Gaussian to integrated line emission in (m, v) plane. The typical statistical error is 

0.01 Jy km s"', but the calibration uncertainty is roughly 20%. 

'*' Offset from position of continuum fitted elliptical Gaussian. Used in Fig.[3]to plot the positions of the lines. 

(<:) FWHM from Gaussian fit in the {u, v) plane. 

(rf.c) Blended pair, values given for the combined integrated emission. 

*^' Upper limit estimate. 
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4.2. Water associated with tiie outflow in IRAS2A 



0.2 0.0 -0.2-0.4-0.6-0.8 

RA Offset (") 

Fig. 3. Positions of the lines from circular Gaussian fits in the 
(Mv)-plane to the integrated line emission. The positions are plot- 
ted relative to the continuum peak position (cross) and the error 
in position for the continuum peak position is roughly the extent 
of the cross while the extent of the ellipses shows the eiTor in 
position for the lines. Ellipses with the same colors represent the 
same molecule, while the numbers refer to a certain molecule 
and transition which has the same number in Fig. [1] Table |3] 
Table lA.fl and Table |A] 2. Many of the detected lines are system- 
atically offset from the continuum peak. 



here are not aligned with any larger scale outflow, and because 
a tentative velocity gra dient perpendicular to the outfl ow is seen 
towards IRAS4B ( J0r gensen & van Dishoeckll2010ah . we argue 
that the lines detected here do not have their origin in the out- 
flow. The line widths of 1.6 - 4.1 km s ' observed in these data 
(Sect. l4.4T i are also signific antly less than the outfl ow widths seen 
in Herschel-HIFI spectra dKristensen et alJ[2010h . 



Whereas the bulk of the Hj 



O emission toward the three sources 
is associated with the central protostellar component, some ex- 
tended emission at the velocity of the Hl^^O line is seen ex- 
tending toward the southwest in the integrated intensity map for 
IRAS2A (Fig.|2]i. Previous observations of IRAS2A have shown 
that it is oriented such th at the outflow lob e extendin g in the SW 
direction is blueshifted (Ij0rgensen et al.l 12004: J0rgensen et al.l 
l2007h . 

To investigate this emission further, continuum subtracted 
channel maps between 2 and 33 km s"', integrated in steps of 
2. 1 km s"' , are shown in Fig. |4] Starting at the higher velocities, 
coiTesponding to SO2 emission in Fig.[T] strong on-source emis- 
sion is seen, but also the extended outflow emission. The gray 
lines trace the rough position of the peak outflow emission in the 
frames 5.0 to 32.0km s', corresponding to a linear velocity gra- 
dient characteristic for outflow affected material: the peak of the 
emission in each frame moves further away along the outflow 
when moving to more blue-shifted velocities. At 19.5 km s"' 
the knot identified as SO2 has moved about 3" away from the 
central source. The knot persists at this position (3~4" away) 
almost throughout the rest of the maps and at a velocity of 5.0- 
7.1 km s"' it weakens. At 9.2 km s"' emission is seen to emerge 
close to the protostar, which is identified as Hj**©. It strengthens 
slightly and moves away from the source with lower velocities. 
Just as with the SO2 emission it can be traced, although over 
a smaller number of channels and spatial extent (red line). The 
outflow water emission is extended and weak, but shows that 
water is present in the outflow of IRAS2A. 

The detection of extended outflow Hl^O and SO2 emission 
is further strengthened when looking at spectra at different posi- 
tions in the outflow. Spectra extracted at two positions (averaged 
over l"xl" regions) from the low surface brightness material 
in the outflow are shown in Fig. |5] In Fig. |2]the positions are 
marked with a (magenta cross) and b (green cross). Outflow po- 
sition a shows one SO2 component with a weak blue tail and a 
broad Hl^^O component that is blueshifted 4-6 km s ' from the 
rest velocity. This blueshift coincides with the velocity of the 
CH3OCH3 lines 3 and 4, but the absence of the CH3OCH3 fines 
1 and 2, which are detected in the on-source spectra, shows that it 
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Fig. 4. Continuum subtracted channel map between 2 and 
33 km s ' for IRAS2A, velocity in km s ' given in each 
frame, channel width is 2.1 km s '. Contours in steps of 3cr 
(24 mJy beam"' km s"') from a 3cr level. The beam size is given 
in the lower left comer of the first frame. Gray and red broad 
lines trace the different outflow components of SO2 and H^^O, 
respectively . 



is indeed H'^O. At position b only SO2 with a broad blueshifted 
component and a narrow component is visible. 

Thus the extended outflow emission seen towards IRAS2A 
in Fig.[2]is at least partly blueshifted Hj^O in the inner part, close 
to the source and SO2 in the outermost part of the outflow. 



4.3. Constraints on geometry and dynamics 



I Watson et al.l (l2007h suggested that the difference between 
IRAS4A and IRAS4B could be attributed to either hot water 
only being detectable during a short phase of the life cycle of 
YSOs or alternatively a particularly favorable line of sight in that 
system. Because of the lower dust opacity, protostellar envelopes 
are optically thin at submillimeter wavelengths, which allows 
detections of the water emission lines from the smallest scales 
towards the center of the protostars. The detections towards three 
of three targeted sources (|j0rgensen & van Dishoeckl2010al and 
this paper) therefore argue against the first of these suggestions. 
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Fig. 5. IRAS2A off-position spectra, averaged over l"xl" re- 
gions around the positions marked with a and b in Fig.|2]in the 
blueshifted outflow from the imaged data. The spectra have been 
rebinned to twice the velocity resolution (width of 0.23 km s"'). 
The systemic velocity of the Hj^O line is drawn with a vertical 
dotted line at visi - 7 km s"', and the gray outline in the back- 
ground is the on-source spectrum averaged in a region of the 
same size as the other spectra. 



The detected velocity gradient in IRAS4B - which is ab- 
sent in the other sources (Fig. |2]i - may be a geometrical effect. 
A simple solution could be that the sources without a velocity 
gradient are observed close to face-on. This would require that 
IRAS4B with its detected velocity gradient is observed edge-on, 
however. This is not su pported by observations of the outflows 
(|j0rgensen et al.ll2007h : all three sources show well-separated 
red and blue outflow cones suggesting at least some inclination 
of the outflows relative to the line-of-sight. Also, of the three 
sources, IRAS4B in fact shows the most compact outflow: if one 
assumes similar extents of the outflows this would suggest that 
IRAS4B rather than the other sources are seen more close to 
face-on. 

Another way to deduce the inclination of the outflow with 
respect to the plane of the sky is through observations of maser 
spots in the outflows of the sources, although they do not seem 
to agr ee with obser v ations of the larger scale outflow in many 
cases. iMarvel et all ( |2008|) conclude that in IRAS4A-NW, the 
maser emission cannot be directly associated with the large scale 
outflow. For IRAS 4B the masers have been assoc iated with the 
outflow (iMarvel et al. 2008; Des murs et al.ll2009l) . but the posi- 
tion angle of the "maser outflow" still differs from that of the CO 
outflow by some 20-30°. 

We propose that the compact Hj^^O emission has its origin in 
a flattened disk-like structure where inward motions is still dom- 
inating the dynamics, and where the viewing angle defines the 
observational characteristics (Fig.|6]for a schematic drawing). In 
fact, Brinch et al. (2009) suggest, based on modeling of submil- 
limeter observations, that IRAS2A indeed has such a structure 
on a few 100 AU scales. Observing a source with this structure at 
different viewing angles would potentially produce the observed 
differences. 

In the proposed solution, IRAS4B is viewed at an angle, 
closer to face-on, causing the mid-IR emission observed by 
Spitzer to escape from the envelope through the cavities pro- 



M. V. Persson et al.: Subarcsecond resolution observations of warm water towards three deeply embedded low-mass protostars 






Fig. 6. Schematic drawing of a disk dominated by inward mo- 
tions rather than rotation in different observation geometries. In 
ca.se A (IRAS4B), the observed line is narrow and a small veloc- 
ity gradient is seen. In ca.se B (IRAS2A and IRAS4A-NW), the 
observed line is broader and no velocity gradient is observed. 



duced by the outflow, where the optical depth alo ng the line of 
sight i s lower in agreement with the suggestion bv lWatson et al.l 
( I2OO7I) (case A in Fig. |6]l. In IRAS2A and IRAS4A-NW, the 
viewing angle is closer to edge-on, causing mid-IR emission to 
be blocked by the optically thick envelope (case B in Fig.|6ll. 

This solution could also explain the velocity gradient ob- 
served towards IRAS4B and the broader line widths towards 
IRAS2A and IRAS4A-NW. In disks dominated by inward mo- 
tions rather than rotation, the systems viewed at an angle exhibit 
narrow lines and a velocity gradient (case A in Fig.|6), while sys- 
tems viewed edge on show broader lines and no velocity gradi- 
ent (case B in Fig.|6]l. The dynamical structure of the disk would 
move towards a rotationally supported disk on a local dynam- 
ical time-scale, wh ich is a few xlO^ yr at a radius of 110 AU 
dBrinch et al.ll2009l) . 



4.4. Non-masing origin of H\^0 emission 

We now turn the attention to chemical implications of the pre- 
sented data: the subsequent analysis of the data focuses on the 
abundance and excitation of water and other species detected in 
the spectra. 

The primary concern about the Hj^^O analysis in this context 
is whether the Hj^O line could be masing. Many lines of wa- 
ter observable from the ground are seen to be masing, as is the 
case for the same transition stud ied here (3|^ - 2 ?n) but in the 
main isotopologue, H^^O tCernicharo et al.l [l 9941) . Several ob- 
servables argue against that the H2^0 3 1,3 - 22,0 hne is masing in 
the sources investigated here. 

First, the line widths observed in the central beam to- 
wards the continuum peak are comparable for all the detected 
species in each of the sources, as was also seen in IRAS4B 
(J0rgensen & van Dishoeck 2010a b). Although the mean line 
width is somewhat broader in IRAS2A and IRAS4A-NW com- 
pared to IRAS4B they are still significantly more narrow towards 
the central source (<5 km s-') than in the broad outflow com- 
ponents seen in the Herschel data (> 5 km s-') towards these 
sources (Kristensen et al. 2010). For IRAS2A, the mean FWHM 
and standard deviation of the fitted lines is 2.8+0.7 km s-', for 
IRAS4A-NW 2.4+1.0 km s-i, and for IRAS4B l.l±0.2km s-'. 



In addition to this, the variations in intensities between the dif- 
ferent lines in each source are smaller than the variations in line 
intensity between the sources. 

In IRAS2A, iFuruva et alj ( l2003h report that maser emission 
at 22 GHz is highly variable and has a velo city of 7.2 km s' ', i.e. 
marginaUy red-shif ted. For IRAS4A, Fur uva et aH (l2003l) and 
iMarvel et al.l (l2008l) show that the strongest maser component 
lies at ~I0 km s~' where we do not detect any strong H^^O mil- 
limeter emission. 

Taken together, these arguments suggest that the H^^O emis- 
sion detected towards the central positions of the sources is not 
masing, but instead is thermal emission with the same origin as 
the other detected species. 

4.5. Column densities and abundances determinations 

Having established that the Hj^O line is unlikely to be masing, 
the column densities for each species can be determined and 
compared directly. The emission from the different molecules 
is assumed to be optically thin and to uniformly fill the beam 
(valid since the deconvolved size in Table[3]is comparable to the 
beam size). The excitation is taken to be characterized by a sin- 
gle excitation temperature Tex, which does not have to be equal 
to the kinetic temperature. 



4.5.1 . Water column density and mass 

Table|4]summarizes the inferred column densities of the various 
molecules assuming rex=170 K. The gas-phase H2O abundance 
and mass are calculated with an isotopic abundan ce ratio of 
i6q^i8q of 560 in the local interstellar medium (Wils on &Roodl 
1994). The masses of gas phase H2O in the sources are given 
in the bottom of Table |4] in units of Solar and Earth masses, 
along with the fraction of gas phase water relati ve to H2. As 
pointed out bv'j0rgensen & van Dishoeckl (l20I0ah . in IRAS4B, 
the derived column density of H?,0 is almost 25 times higher 
than deduced bv lWatson et al] (l2007h . The deduced mass of the 
gas phase water in IRAS2A is roughly six times higher than in 
IRAS4A-NW and 90 times higher than in IRAS4B. Thus, wa- 
ter is present in the inner parts of all of the sources and the 
non-detection in IRAS4A and IRAS2A at mid-infrared wave- 
lengths is due to the high opacity of the envelope combined with 
a strong outflow con tribution to H^^O lines in IRAS4B (see also 
lHerczegetaI.ll20Ill) . 



Table 4. Derived properties of the detected species over the de- 
convolved sizes of the lines for all sources in the project. 





1RAS2A 


IRAS4A-NW 


1RAS4B 


■'VcHaOCHj" 

A'so," 


1.2x10"" 
2.9xl0i« 
<1.4xl0"' 
1.8x10'* 


6.8x10'" 
7.9x10" 
1.1x10" 
8.2x10" 


5.6x10'" 
3.9x10'^ 
2.8x10"* 
1.3x10'^ 


Mh,o (Mq) 
Mh,o (Me) 


1.8x10-" 

0.61 

4.2x10-" 


2.7x10-' 

0.09 

1.4x10-' 


1.5xlO-« 

0.006 

7.8x10-' 



*"' In cm-- assuming rex=170 K. The typical error is a combination of 
the assumption of the excitation temperature and the flux calibration 
uncertainty (~20%). The determined column density varies with T^x in 
the interval 100-200 K differently for each molecule, see Fig.|8]and the 
discussion in Sect. 14.5.31 *' Using only lines 1, 2, 3 and 4. 
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Fig. 7. Left: Fraction of gas phase water relative t o H2 (from the 
continuum observations of J0rgensen et al.ll2009 ^ for the differ- 
ent sources, including an upper limit estimate for IRAS4A-SE. 
Right: Ratio of the column density of the molecules relative to 
water. The C2H5CN value for IRAS2A is an upper limit. 



The ratio of H2O gas phase abundance to th e total H2 abun- 
dance (derived from continuum observations dj0rgensen et al.l 
I2OO9I) ') in Fig. |7] shows that the fractional abundance of gas 
phase water in IRAS2A is much higher than in IRAS4A-NW 
and IRAS4B. Still, the gas-phase water likely accounts for only 
a small fraction of the total amount of water (gas+ice) on 
those scales. Typical H2O to H2 abundance ratios of 10"^ - 
10'" ^ have been deduced in shocked regions close to protostars 
jHarwit et al. 1998; Nisini et al. 1999 ), and in ices in the col d 
outer parts of protostellar envelopes (iPontoppidan et all |2 0041) . 
In shocks, high temperature gas-phase reactions drive much of 
the oxygen into water (Draine et al. 1983). In more quiescent re- 
gions the ice is expected to evaporate off the grains when heated 
to tem peratures above 10 K in close proximity to the central 
source dFraser et al.ll200lh . The fact that our inferred abundance 
ratios are generally much lower than 10""* suggests that a large 
fraction of the dust mass in the inner ~ 100 AU has a temperature 
below 100 K, with the fraction varying from source to source. In 
colder regions, most of the water gas originates from UV pho- 
todesorption of ices. Assuming a total H2O to H2 abundance ra- 
tio of 10""* (gas-i-ice), the fraction of water in the gas phase is 4% 
for IRAS2A, 0.1% for IRAS4A-NW and 0.01% for IRAS4B of 
the total water content (gas+ice). This difference could be due, 
for example, to varying intensities of the UV radiation from the 
newly formed stars and geometry causing differences in the rel- 
ative amounts of material heated above the temperature where 
water ice mantles are evap orated. Towards th e more evolved 
Class II source TW Hvdra. iHogerheiide et al.l (1201 1) deduce a 
gas phase water to total hydrogen abundance of 2.4x10"** using 
Herschel. This value is of the same order of magnitude as our 
results, and places TW Hydra just above IRAS4B in Fig.|7] 

4.5.2. Other molecules 

In the right plot of Fig. |7] the abundances of the detected 
molecules relative to H2O are presented. The CH3OCH3 and 
SO2 content relative to Hj^O is roughly comparable between 



the different sources. This suggests that for CH3OCH3, SO2 and 
H^**0 (and thus H2O) similar chemical processes are at work in 
the three sources. C2H5CN is the least abundant of the molecules 
for all the sources. In IRAS2A the high abundance of water 
is contrasted by the relatively low abundance of C2H5CN, for 
which only an upper limit is deduced. This relatively low abun- 
dance in IRAS2A is robust even when much lower excitation 
temperatures are assumed (see Sect. I4.5.7t . 

It is also evident from Fig. |7] and Table |4] that the 
CH3OCH3/C2H5CN column density ratio varies significantly 
between the sources. This ratio is 10^ for IRAS2A, 60 for 
IRAS4A-NW and 200 for IRAS4B assuming Tex =170 K. The 
differences in relative abundances of N and O-bearing complex 
organics have be en interpreted as an indicato r for different phys- 
ical parameters (jRpdgers & Chamlev"2001') and differen t initial 
grain s urface compo sition (Charnlev et al. 1992; Bottinel li et al.l 
l2004t) . ICaseIIi et ^.1 (Il993h argue that the (dust) temperature his- 
tory of the cloud, in particular whether it was cold enough to 
freeze out CO, is an important factor for the resulting ratio. 

4.5.3. Effect of excitation temperature 



The results presented in Sections [4.5. H and 14. 5. 2l assume an ex- 
citation temperature of 170 K for all species. In this section we 
demonstrate that all our conclusions are robust over a wide range 
of assumed temperatures. The adopted rex= 170 K is the same as 
found from the model fit to the mid-infrared water lines observed 
by Spitzer ('Watson et al."2007') and allows for easy comparison 
to this paper and J0rgensen & van Dishoeck (2010a). 

Although the origin of the mid-infrared and millimeter lines 
are likely not the same, this temperature is a reasonable first 
guess appropriate of the conditions in the inner 150 AU of proto- 
stellar systems. From modeling of dust continuum observations, 
the density in the inner envelope/disk is found to be high, at least 
10^ cm"^, and the dust temperatur e is thought to be at least 90 K 
for a fraction of the material (e.g., |j0rgensen et al.ll2002[ l2004l) . 
This is also the temperature at which water is liberated from the 
grains; any gas-phase water detected in the high excitation line 
observed here (£'u~200 K) likely originates from the warmer re- 
gions (Sect. l4.5.lT l. 

The current observations and literature data do not pro- 
vide enough information on the line emission at these scales 
to constrain the excitation temperature directly. For none of 
the detected species are there more than two previous mea- 
surements and those detections either have large uncertainties 
due to their weakness ('j0rgense n et al.ll2005h or are sensitive to 
larger scale emission including ou tflow and envelope because 
of the large beam (BottinelH et ak 2004; Botti nelli eFaP 120071: 
Wakelam et al. 2005) . F or the deeply embedde d low-mass pro- 
tostar IRAS 1 6293-2422. lBisschop et all (l2008l) derive rotational 
temperatures between 200^00 K from several compact emis- 
sion lines of complex molecules using SMA data. Furthermore, 
for the outflow toward IRAS4B iHerczeg et al.1 (1201 Ih deduce 
Tex^l 10-220 K with recent Herschel R\CS data. 

Figure [8] explores what happens with the column densities 
at significantly different excitation temperatures than those as- 
sumed here. It shows two different plots for the three sources. 
The plots in the top row show the ratio between the column den- 
sity of the various molecules and that of water, whereas the plots 
in the bottom row show the column densities themselves, all 
as functions of excitation temperature. The ratios (Nx/Nuio) do 
not change significantly between 100 and 275 K. The gray area 
in the lower plots shows the range of excitation temperatures 
for which the column density of water is the same in the three 
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Fig, 8. Plots showing (upper) the ratio between the detected 
molecules and water (H2O), assuming that the species have the 
same Tex and (lower) the column density of all the detected lines 
(including Hj^O). The colors and line style indicate the different 
molecules, the gray vertical line mark T^x- 170 K and the gray 
area in each plot indicate the temperatures where the column 
density of Hi^O is the same for the sources. 



(exp(-£'u/k;rex)) does not compensate for the fall that the par- 
tition function (proportional to Tex*, where b>0) introduces. In 
the temperature range 40 to 250 K the CH3OCH3/C2H5CN col- 
umn density ratio varies with a factor of about 30. Looking only 
at the CH3OCH3/C2H5CN column density ratio, the excitation 
temperature for IRAS2A has to be rex=45 K, and for IRAS4B 
rex=85 K, to match that of IRAS4A-NW. 
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Fig. 9. Plot of the ratio between CH3OCH3 and C2H5CN in 
IRAS2A (left y-axis, filled line), and the column densities de- 
rived for those molecules (right y-axis, dashed/dash dotted line) 
for temperatures between 25 and 275 K. The vertical lines show 
r=40 and 170 K. 



As argued above based on the physical characteristics of the 
sources, values of Tgx significantly below 100 K and large dif- 
ferences in Tex between molecules and sources are not expected 
on the scales probed here. Also, the agreement in widths, extent 
and positions between the lines argues against this. If the temper- 
ature is allowed to vary in the range 100-250 K the conclusions 
drawn in this paper are unchanged. 



sources. Only if the excitation temperature in IRAS4A-NW and 
IRAS4B is lower than 50 K can all sources have comparable 
water column densities based on the observed intensities. 

The different molecules could of course each have differ- 
ent excitation temperatures. Specifically, can different excita- 
tion temperatures be the cause of the large range in inferred 
CH3OCH3/C2H5CN column density ratios? Figure |9]shows this 
ratio and the individual column densities for excitation temper- 
atures from 25 to 275 K for IRAS2A. The corresponding plots 
for the other two sources show similar trends (see also Fig. |8]). 
The temperature does not affect the column density estimate for 
C2H5CN significantly, showing that our upper limit estimate is 
rather robust, and that the differences in this ratio between the 
sources cannot be attributed to differences in the excitation of 
C2H5CN alone. On the other hand, changing the excitation tem- 
perature for CH3OCH3 affects the column density and the ratio 
to a greater extent. This difference in behavior with Tex for the 
two molecules is caused by the transitions' upper energy lev- 
els (£u); 18 K for CH3OCH3 versus 122 K for C2H5CN. For 
CH3OCH3, its lower value of E^, causes the column density es- 
timate to decline more rapidly with Tex than for C2H5CN: the 
contribution from the exponent in the Boltzmann distribution 



4.6. Formation of molecules In the protostellar environment 

The simultaneous observations of water and the complex mole- 
cules can also be used to shed new light onto where and how the 
formation of various molecules takes place in the protostellar 
environment. The basic scenario is that many of the complex or- 
ganic molecules are formed as a result of grain-surface reactions 
during the collaps e of the protostellar envelope s as material is 
he ated slowly ( see Herbst & van Dishoeckll2009l for a review). 

iGarrod & Herbst (2006) use a gas-grain network of reac- 
tions to investigate the formation of certain complex organic 
molecules in hot cores and corinos during the warm-up and evap- 
orative stages. The results show that the chemistry occurs both 
on grain surfaces and in the gas phase, with the two phases being 
strongly coupled with each other. Large abundances of complex 
organic molecules are present at the end of this warm-up stage, 
and show that a period of 10** - 10^ years is vital for the forma- 
tion of complex molecules. 

A key question is thus whether the material in these low- 
mass protostars spends enough time at moderate temperatures 
of 20^0 K to form complex organic molecules during the col- 
lapse of the cloud and formation of the disk. A parcel of ma- 



M. V. Persson et al.: Subarcsecond resolution observations of warm water towards three deeply embedded low-mass protostars 



terial originating in the cold core will pass through the warm 
inner envelope during the collapse and be incorporated into the 
disk before being accreted onto the star. Depending on its precise 
origin, the material can follow many different routes, but explo- 
ration of a variety of collapse models shows that most parcels 
indeed spend the required period at moderate temperatures to 
form the complex organics, either in the inner envelope or in 
the disk (Visser et al. 2009 ). Some of the routes result in the 
material being dir ectly exposed to U V radiation near the out- 
flow cavity walls dVisser et al.l 120121) . After evaporation, gas- 
phase temperatures in t he hot core are generally too low ( i.e., 
T <230K,R> 20 AU. lLarsonll9"7llJ0rgensenet al.l2002l) for 
any significant gas-phase formation of water to take place. The 
fact that we observe narrow emission lines of complex molecules 
on small scales confirms that some fraction of the collapsing 
cloud has spent enough time under conditions that lead to the 
formation of complex species that will enter the disk. 



5. Summary and Outlook 

We have presented subarcsecond resolution observations of the 
Hj^O 3 13 -22,0 transition at 203.4 GHz towards a small selection 
of deeply embedded protostars in the NGC 1333 region. In ad- 
dition to the targeted H^^O molecule, SO2 and complex organic 
molecules are seen. 

- Water is successfully detected towards IRAS2A and one 
component of the IRAS4A binary. These data add to the pre- 
vious detect ion of the 3 1.3 - 22.0 transition of H J,^0 towards 
IRAS4B bv lJ0rgensen & van Dishoeckl (l2010al) . The origin 
of the emission differs from that seen with both Herschel and 
Spitzer, where larger-scale emission (protostellar envelope, 
outflow shocks/cavity walls) is detected. 

- The gas-phase water abundance is significantly lower than 
a canonical abundance of lO"'*, suggesting that most wa- 
ter is frozen out on dust grains at these scales. The frac- 
tion of water in the gas phase is 4% for IRAS2A, 0.1% for 
IRAS4A-NW and 0.01 % for IRAS4B of the total water con- 
tent (gas-Hice). 

- IRAS2A shows H^^O eission not only at the source position, 
but also in the blueshifted outflow. In IRAS4A we only de- 
tect water towards IRAS4A-NW. None of the lines presented 
for the other sources are detected towards IRAS4A-SE. 

- Based on the small species-to-species variations, in both line 
strength and line width, we argue that the Hj^O 3 13 - 22,0 
emission line is not masing in these environments and on 
these scales. 

- While IRAS4B shows narrow e mission lines and a tentative 
velocity gradient, reported by lj0rgensen & van Dishoeckl 
(l2010al) . the other sources show broader lines and no veloc- 
ity gradients. Additionally, IRAS4B is the only source of the 
three detected in mid-infrared water lines by IWatson et al.l 
(12007.) . We propose that a disk-like structure dominated by 
inward motions rather than rotation with different orienta- 
tions with respect to the observer could explain the observed 
line characteristics. 

- The abundances of the detected molecules differ between 
the sources: the abundance of C2H5CN with respect to wa- 
ter, compared to both SO2 and CH3OCH3, is much lower 
in IRAS2A than the other sources. This shows that the for- 
mation mechanism for C2H5CN is different from the other 
molecules, indicating that it is not trapped in the water ice 
mixture or related to the sublimation of water 



The detections of Hj^O emission from the ground show 
good promise for using Hj^O as a probe of the inner regions 
of protostars and for studying the formation of complex organic 
molecules in the water-rich grain ice mantles. Higher resolution 
observations, using the PdBI in the most extended configura- 
tion, or with the Atacama Large Millimeter/submillimeter Array 
(ALMA) band 5 could be used to reveal the kinematical struc- 
ture of the 10-100 AU scales and the dynamics of the forming 
disks. Observations of higher excited lines of Hj^O and the com- 
plex organics in other ALMA bands (e.g., band 8) could provide 
a better handle on the excitation conditions of water on few hun- 
dred AU scales in the centers of protostellar envelopes. Finally, 
observing more protostellar sources in the Hl^^O 3 13 - 22,0 tran- 
sition would provide a better statistical sample for comparisons 
and conclusions. 
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Appendix A: Tables 



Table A.l. Tabulated data for the different lines 
detected frorn JPL (Pickett et al. 1998) and CDMS 
(lMulleretal.ll200l . 

No. Molecule Transition RestFreq. £„ logA,^ 

(GHz) (K) (log.r') 



1 


CH,OCHi 


33 ~ 22 1 EE 


203.4203 


18.1 


-4.61 


2 


CH,OCHi 


330-22 1 AA 


203.4187 


18.1 


-4.23 


3 


CH3OCH3 


33.0-22.1 AE 


203.4114 


18.1 


-4.23 


4 


CH3OCH3 


33.1 -22.1 EE 


203.4101 


18.1 


-4.46 


5 


Hi«0 


3 1,3 -22,0 


203.4075 


203.7 


-5.32 


6 


CH3OCH3 


33.1-22., EA 


203.4028 


18.1 


-4.39 


7 


C2H,CN 


232,22-222,2, 


203.3966 


122.3 


-3.15 


8 


SO2 


12o,12 - 11|,1, 


203.3916 


70.1 


-4.07 
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Table A.2. Position and FWHM of the detected lines from Gaussian fits to the spectra. 









1RAS2A 






1RAS4A-NW 




No. 


Molecule 


Position 


Width 


Interval" 


Position 


Width 


Interval" 






(kms-') 


(kms-') 


(kms-') 


(km s-') 


(kms-') 


(kms-') 


1 


CH^OCH, 


-11.7±0.1 


2.3 ±0.1 


-14.0- -7.0 


-11.3±0.1 


1.6 ±0.1 


-13.0 --7.3 


2 


CH,OCHi 


-9.3 + 0.1 


2.3 ±0.1 


- 


-9.0 ±0.1 


1.6 ±0.1 


- 


3 


CH^OCHi 


+ 1.3 ±0.1 


2.3 ±0.1 


-1.0 -+4.4 


+ 1.7±0.1 


1.6 ±0.1 


0.0 - +4.7 


4 


CH^OCH, 


+3.2 + 0.1 


2.3 ±0.1 


- 


+3.6 ±0.1 


1.6 ±0.1 


- 


5 


m^'o 


+7.4 + 0.1 


4.0 ±0.1 


+4.4- +11.4 


+6.6 ±0.1 


2.9 ±0.3 


+4.7 - +9.6 


6 


CH,OCH, 


+ 15.2 ±0.2 


2.3 ±0.1 


+ 12.9- +17.5 


+ 14.4 ±0.1 


1.6 ±0.1 


+ 12.0- +15.2 


7 


CjH^CN 


- 


- 


- 


+22.8 ± 0.2 


3.8 ±0.4 


+ 19.0 -+26.9 


8 


SO, 


+30.4 ±0.1 


3.8 ±0.1 


+26.6 - +33.9 


+30.4 ±0.1 


4.1 ±0.2 


+26.9 - +34.6 



'"' Interval over which the integrated line intensity (Table |3) is estimated. 
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